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Abstract. Recent observations suggest molecular line ratios in millimeter and 
submillimeter bands may be a good tool to reveal the long-standing question on 
the origin of energy sources in obscured active galaxies - AGN and/or starburst. 
Observations of actual molecular medium show in general inhomogeneous struc- 
tures as well as high-resolution hydrodynamic simulations do. In order for pre- 
cise interpretation of emergent line emission from the inhomogeneous molecular 
gas to probe the dominant energy source of active galaxies, we study character- 
istic features of emergent intensities via three-dimensional non-LTE (non-local 
thermodynamic equilibrium) line transfer simulations. Our results succeeded 
in making clear 1) the necessary conditions for HCN/HCO + -dichotomy, and 2) 
importance of dumpiness on intensity ratio and its interpretation. These results 
are obtained for the first time by our realistic three-dimensional simulations, 
and line transfer simulations will be a powerful tool to comprehensive studies of 
extragalactic interstellar medium (ISM) in forthcoming ALMA (Atacama Large 
Millimeter/submillimeter Array) era. 



Introduction 



Molecular gas in external galaxies is a subject of crucial importance for obser- 
vational and theoretical studies of galaxy formation. Among the interstellar 
medium (ISM), compact molecular gas around an active galactic nuclei (AGN) 
is expected to be energy budget from AGN and/or the nuclear starburst of 
possible relevance. Recent observations suggest line ratios in millimeter and 
submillimeter bands may be a good tool to reveal the long-standing question of 
the origi n of activities in obscured active galaxies - AGN or nuclear starburst 
(see e.g. iKohno et al]|200ll : iKohnol [20051 ; llmanishi et~alll2007l ). 

Current observations present compact molecular gas of size ~ 1 kpc or less, 
though, hydrodynamic simulations predict highly inhomogeneous and clumpy 
substructures in a " compact" molecular gas at the centers of active galaxies 
( Wada &: Normanl l200ll ; IWada <fe Tomisakal 12005). In such an inhomogeneous 
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molecular gas, energy level populations of emitting particles will be also com- 
plex and complicated one, which makes derivation of correct physical features of 
emitting ISM from emergent line intensities difficult. Then we study excitation 
conditions and line intensities from an inhomogeneous molecular gas at the cen- 
ter of an active galaxy, with three-dimensional non-local thermodynamic equilib- 
rium (non-LTE) line trans fer simulations along with the results of high-resolution 
hydrodynamic simulation ( Wada Sz Tomisa ka 2005 ; Yamada. Wada Sz Tomisakal 
20071 : 1 Yamada & Tomisakal I2008D . 



2. Numerical Calculations 

We first performed hydrodynamic simulations of a compact (R < 32 pc) ISM in a 
steady gravitational potential of supermassive blackhole (Msmbh = 1O 8 M ) and 
host galaxy. And then we performed non-LTE line transfer simulations as a post 
process using a snapshot result of hydrodynamic simulation. We solved non-LTE 
rate equations in statistical equilibrium up to J = 10 and integrated a standard 
radiative transfer equation along all the sampling rays to calculate mean in- 
tensity J = l/(47r) / I u dQ. We calculated these two equations iteratively until 
energy lev el distributions and mean intensity field c onverge with ~ 10~ 6 relative 
preci sion (jHogerheijdge van der Tald [2000 ; see lYamada. Wada k, Tomisakal 



20071 for details). 



3. HCN and HCO+ Dichotomy of Active Galaxies 



Since the pioneering work of lKohno et al. ( 200ll ). it has been argued that inten 



sity ratio of HCN and HCO + rotational lines would be able to probe obscured 
cen tral energy sources of active galaxies in terms of ch emical abundances (see 
e.g. lKohnoir2005l ; llmanishi et al. 2007 ; Baan et al. 20081 ) . The basic idea under- 



lying this dichotomy is that strong X-ray from AGN will form X-ray dominated 
regions (XDR) compared with stellar feedback from nuclear starburst that will 
form photodissociation region (PDR), and different chemical structures in XDR 
and PDR imprinted in molecular line ratio indirectly enables unveil the central 
energy sources. 

In Figure Q] we show our numerical simulation results with an assumption 
of spatially uniform chemistry. Distribution of integrated intensity of HCN 
(1-0) line (Fig. QJci)) becomes quite inhomogeneous, reflecting the inhomoge- 
neous structure in the molecular gas model. Numbers of bright regions of which 
radius is ~ lOpc or less are observed. Current millimeter observational instru- 
ments cannot resolve these small structures except for the closest galaxies, but 
forthcoming ALMA (Atacama Large Millimeter /submillimeter Array) telescope 
will reveal these internal structures in a "compact" gas at the center of a dis- 
tant galaxy. In order to examine intensity ratio of HCN(l-O) and HCO + (1-0), 
we computed HCN and HCO + lines for a wide range of molecular abundances 
(10~ n < y < 10~ 7 , where y is relative molecular abundance to H2). In Figure Q] 
(b) we plot the ratio of resultant intensities averaged over the " field-of-view" of 
the simulation as a function of y for face-on view. Our results show that for the 
same value of y, HCO + (1-0) is always brighter than HCN (1-0), so that a high 
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Figure 1. Left panel : integrated intensity of HCN (1-0) line in a face-on 
view, right panel : average integrated intensity (J) for HCN (black line with 
crosses) and HCO + (blue line with squares). 



intensity ratio -Rhcn/hco+ greater than unity requires overabundance of HCN 
compared with HCO + (yucN ~ 10 x yHCO+) on average. 

Current models of XDR chemistry predicted under-a bundance of HCN 
(|Meiierink fc Spaansl [20051 : iMeiierink. Israel fc Spaansll2007l b and then original 
picture of HCN/HCO + diagnostics in terms of XDR/PDR chemistry would need 
revision and sophistication. Extended observations have given rise discussion 
on alternative mechanisms that determine -Rhcn/hco+ such as mid-IR photon 
pumping ( Aalto et al. 20071 ). Extended line transfer simulations including non- 
uniform chemical structure and/or continuum pumping will serve as a powerful 
tool to more detailed studies as well as more sophisticated chemistry models 
with realistic gas geometry. 



4. Intensity Ratio of Millimeter and Submillimeter Lines 

Recent progress of submillimeter observations of high- J molecular lines opened 
a way to examine physical conditions of emitting ISM in more detail. Since crit- 
ical densities for LTE population and wavelengths are dependent on J, average 
intensity ratio will be affected by density and temperature structures, and the 
effect of inhomogeneity is expected to be stronger for high density tracer because 
of small fraction of dense and thermalized gas. We examined the ratio of J =4- 
3 and 1-0 lines of HCN molecule as representative dense gas tracers with lin e 
transfer simulations described in §2 ( Yamada 2008; Y amada Tomisaka 2008). 



Figure (a) shows averaged intensity ratio as a function of viewing angle 9 
(8 = 0° is face-on view, and 90° is edge-on view). It shows that average #43/10 
decreases with increase of y, which contradicts one-zone estimation that -R43/10 
approaches unity as y (or optical thickness) increases. This result is likely to 
come from multi-phase nature of emitting ISM. 
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Figure 2. Panel (a) : average intensity ratios -R43/10 of HCN for four abun- 
dance values taken from line transfer simulations. Panel (b) : comparison 
of analytic formula (solid line) for a simple isothermal two-phase model and 
numerical line transfer results of isothermal ISM models. In panel (c) numer- 
ically evaluated -R43/10 for an isothermal two-phase modelling are plotted as 
a function of r\ and Tki n , and a fixed average density n avc = 100 cm~ 3 . 



4.1. Two-Phase Modelling and Analytic Formula of -R43/10 

We model multi-phase ISM as a simple two component one based on the knowl- 
edge of excitation mechanisms. Since we know that average kinetic temperature 
of model ISM (174 K) is sufficiently higher than the excitation energy of J = 4 
level of HCN, we assume single temperature ISM. Furthermore in order for an- 
alytic examination we assume optically thin over a whole region. Then our 
two-phase model consists of two density components, one is dense (n > n cr it) 
clumps of which volume fraction is 77, and another is tenuous (n <C re C rit) ambi- 
ent. 

Considering that in dense (n > n CT - lt ) regions level population is close 
to LTE, and in warm {T\^ m 3> /^j,j-i/&b) and tenuous (re <C re C rit) regions 
level population is a balance of spontaneous decay and collisional excitation 
(|Yamada et al.l 12007). average intensity ratio -R43/10 is described as a ratio of 



sum of intensities from both regions 

R _ Evol(^H 2 j//4^43^43)^43 2 _ ( x Svol (^H 2 /4^43 ) 

43/10 Evol^Hay/l^lO^lO^lo 2 V^IO/ Evol( n H 2 /l^lo) 



K 7743QH 2 /4^43)d + i 1 ~ ^3) Qh 2 f\Mz)t ^ 

Vw( n H 2 fiA w )d + (1 - ?7io)(«H 2 /i^io)t' 

Figure [2](c) shows an numerical evaluation of thus obtained -R43/10 It is obvious 
that average -R43/10 is quite sensitive to n, and can take almost any value up 
to theoretical maximum (=16 in this case) even if average density (n ave = 100 
cm~ 3 ) is smaller than critical density (re cr it = 10 5 cm -3 for J = 1-0 line). In 
other words, intensity ratio close to unity does not guarantee a thermalized 
population. 

In Figure[2](b) we plot numerical evaluation of equation (pTJ) with a definition 

Of r/43 = 7710 = ??sim 

= jVgrids(^i/ra > egj/goexp(-^i /fc g T kin )) 

' grids 
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Compared with results of line transfer simulations of isothermal hydrodynamic 
data generated by original input model, we can see that the analytic formula 
of -R43/10 (H|) can fit numerical simulation results within a factor of two in spite 
of simplicity of this formula. We thus concluded that multi-phase modelling 
based on excitation condition considerations along with numerical line trans- 
fer simulations could help construct more sophisticated ways to derive correct 
physical features and structures from observed line intensities. Our results also 
demonstrated the importance of spatial structure of density (and temperature) 
on emergent intensity besides chemical structure. 

5. Summary 

We showed theoretical considerations on molecular line diagnostics of active 
galaxies by numerical experiments of three-dimensional line transfer simula- 
tions. Two kinds of results (-^hcn/hco -1- an< ^ ^43/10 of HCN) showed that simple 
models need revisions to take into account of inhomogeneities of emitting ISM. 
Numerical experiments can achieve a high-resolution examination in prior to 
the forthcoming ALMA project, and will be a powerful tool to examine physics 
imprinted in emergent intensities. 
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